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The quasi-classical model in a gauge theory with the Yang-Mills (YM) field is developed. On a 
basis of the exact solution of the Dirac equation in the SU{N) gauge field, which is in the eikonal 
approximation, the Yang-Mills (YM) equations containing the external fermion current are solved. 
The derived solutions are quantized in the quasi-classical approach. The developed model proves 
to have the self-consistent solutions of the Dirac and Yang-Mills equations at A'^ > 3. Thereat 
the solutions take place provided that the fermion and gauge fields exist simultaneously, so that 
the fermion current completely compensates the current generated by the gauge field due to it 
self- interaction. The obtained solution are considered in the context of QCD. 

PACS numbers: 11.15.-q, 03.70.-|-k, ll.lO.-z 

d 
o 

^ • I. INTRODUCTION 

Q 

CO ' A study of non-Abelian gauge fields plays an important role in the modern field theory The non-Abelian 
gauge field are a basis of QCD jj. The knowledge of solutions of the YM equations enable us to understand features 
of processes in the strong interacting matter generated in collisions of heavy ions of high energies @. Primarily, this 
concerns studying the observable states of such matter as well as the processes accompanying evolution of the medium. 
O I Studying the non-Abelian gauge fields has a very long history which started by the classic paper by C.N.Yang and 
1^ R.L.Mills[6]. Since the paperQ has issued a lot of papers0-[l3 have been devoted to deriving the solutions of the 
. YM equation in various situations. The solution of the source-less YM equation in terms of plane waves was derived 
' in0,Q. A wide class of solutions of the YM equation concerns (1+3) Minkowski space-time in the presence of external 

sources [ol-fTsj. The YM equation were solved for the SU{2) gauge field. The spherical symmetric solutions are 
J> ' found for the SU{2) fields in some specific case of (1+2) space-time in Ref . jisl . The Dirac equation in the presence 
0^ \ of the SU (3) YM field is considered[l^ [l3 terms of the confinement problem. The quark confinement in the curve 
■ space-time is studied in Ref. [l8l| . Rather detailed review of the paper devoting to quantizing the YM field is in the 
monographs by A.Slavnov and L.Faddev[l9j. 

The consistent consideration of the strong interacting particles (generated, for example, in collisions of high energy 
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ions), generally, demands solving the Dirac and Yang-Mills equations simultaneously. The first step in studying such 
problem, naturally (see Ref. [20l. l2l| ) , is an attempt to derive the solution of these equations when the YM field has 

^— ' the form of some modified plane wave so that both the Dirac and Yang-Mills fields will be in the confined region of 
. . ' space. The knowledge of the self-consistent solution of the Yang-Mills and Dirac equations in such approximation 
K>" ■ allows us to obtain the exact Green's function of a fermion field. As a result, it enables to drive both the renormalized 

k>( ' vertex functions and effective mass of a fermion as well as to calculate the observable characteristics of the strong 
interacting matter generated in collisions of high energy ions[5|] beyond the perturbation theory. 

In the present paper the quasi-classical model in the SU {N) gauge theory with the Yang Mills field is developed. 
The self-consistent solutions of both the nonhomogeneous Yang- Mills equation and Dirac equations in an external 
field are derived when the gauge Yang-Mills field is in the eikonal form. The obtained solutions are quantized in the 
quasi-classical approximation. It is shown that the self-consistent solutions of such equations take place when > 3. 
They occur provided that the fermion and gauge fields simultaneously exist, so that a fermion current completely 
compensates the current generated by the gauge field due to it self-interaction. Thereat, there is no energy flux from 
the range of space where the fields are localized. 

In the context of the multi particle problem in matter, the derived solutions mean that the Yang-Mills field is the 
modified circularly-polarized wave which intensity depends strongly on fermion density. In this way, there are both 
the individual and collective states of fermi-particles in matter. When the matter is in equilibrium the type of the 
fermion states depends strongly on such parameters as the temperature and density of the medium. The fermion states 
appear to be one-particle at rather large temperature. With decreasing the matter temperature they are rearranged 
so that the collective states of the fermions interacting with the YM field arise. Thereat, interaction between the 
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fermions and YM field leads to the re-normalization of a fermion mass in the mean. The re-normalized mass depends 
significantly on the temperature of the matter. 

The problem of the hadronization of an equilibrium quark-gluon plasma is considered. In the case of the hadroniza- 
tion into the lightest hadrons due to the phase transition of the first kind a mass of the hadron is calculated. It 
appears to be of the order of the mass of a free pion. 

The paper is organized as follows. Sections II consists of the statement of the problem and the description of the 
main approaches. In Section III the Yang-Mills equation is studied when the gauge field has the form of the eikonal. 
Section IV is devoted to deriving the exact solution of the Dirac equation in the external SU{N) gauge field being in 
the form of the eikonal wave. The quantum consideration of the problem is developed in Section V. The specifics of 
the Yang-Mills and fermion fields obtained in the framework of the developed model are studied in Section VI. In the 
section VII the derived results are considered in the context of QCD. The conclusion is the last section. Appendix A 
contains the calculation of the operator exponent appearing in solving the Dirac equation in an external non-Abelian 
field. The detailed calculations of some important integrals are in Appendix B. 

II. THE YM EQUATIONS IN THE PRESENCE OF EXTERNAL CURRENT 

We consider the SU{N) gauge field A'^ generated by a fermion current. It satisfies the following equations [T9I [2^ : 



d.F^^ix) ~ a ■ f,,^Al{x)Frix) = -gJa-{x) (1) 

F^ix) = d'^A^.ix) - d'^A^ix) ~ g ■ fJ^AUx)A^A^), (2) 

Ja'ix) = *(a;)7''ra*(x), (3) 
where the fermion fields 5* (a;), ^(o;) are governed by the Dirac equation: 

{ind^ + ig-A';^ix)Ta)-m}M'{x)^0 (4) 

^{x){ir{'d^-tg-A*l{x)Ta)+m}=0; = ^A,. (5) 

Here, m is a fermion mass, g is the coupling constant; 7" are the Dirac matrixes, x = x^ = (x'^;x) is a vector 



in the Minkowski space-time; = (5/9t; V); the Roman letters numerate a basis in the space of the associated 
representation of the SU [N) group, so that a,b,c— 1 . . . — 1. We use the signature diag (C/'^") = (1; — 1; — 1; — 1) 
for the metric tensor Q^'^ . The line and "dagger" over ^> mean the Dirac and hermitian conjugation, respectively i23|]. 
Summing over any pair of the repeated indexes is implied. 

The symbols Ta in Eqs.(3)-(5) are the generators of the SU(N) group which satisfy the commutative relations and 
normalization condition: 

[Ta, Tfc]„ = TaTt - TtTa = lU'T,; = -2 z Tr ([T,, n]_ T,) (6) 

Tr [Ta n) = ]^5ab; (7) 

where /^^"^ are the structure constant of the SU (N) group, which are real and anti-symmetrical with respect to the 
transposition in any pair of indexes; 6ab is the Kroneker symbol. In the matrix representation the operators (2 Ta) 
coincide with the Pauli and Gill-Mann matrixes when TV is equal to 2 or 3, respectively. 
It directly follows from Eqs.(6), (7) that 

[Ta, Tfc]+ = Tan + TbTa = ^5ab + dabcT' , dabc = 2 Tr {[Ta,Ti,]^ T,) (8) 

where dabc is real and symmetrical with respect to the transposition in any pair of indexes. 

The main goal is to derive the self-consistent solutions of of Eqs.(l)-(5) which will be localized in the confined region 
of space. We find the solution when the field ^^^(x) is in the form: 
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Al[x) = A'^^i^ix)), (9) 
where '^{x) is some scalar function in the Minkowski space-time which is such that: 

(a^(^)(a^(^) = fc^fc^ = o; (10) 

The last formula determines the well known eikonal approximation where i^ix) can be interpretable as the function 
governing the wave surface of the field A^^. 
We take the axial gauge for the field ^Ji(a;) : 



^^'Al = 0; k^'Al = 0, (11) 



where the dot over the letter means differentiation with respect to the introduced variable Lp. 

Taking into account of both the dependence of A'^{x) on the variable x via the function Lp(x) and formulae (10), 
(11), we derive from Eqs.(l), (2): 

25^'^ fa,'Al{^)Al{^) - (d^dM^)) ■ A: - gk- /,,-4(^)i^f(^) +gVa6^ fc {4(^)Ar(^)A$f (^)} = -g J^ix)- 

J^ix) = ^ifixh'^Ta'ifix). (12) 

It follows from Eq.(12) that in order to derive the solution of the YM equation it is necessary to calculate the 
fermion current Ja{x), which, in its turn, is governed by the solutions of the Dirac equation in the external field 

To do it we assume that the field A„(i^) can expanded as follows in the local frame: 



=^(4i)(¥') cos (^(a;) +¥>«)+ e'(2^{ip)sin{ip{x)+ipa)) + Ba d''ip{x) 
e('i)e(2)^ = e^i)fc,. = 6^2)^1' = 0; e'(^^ = e'(^y, e\^.^=-e\^y k" = d" ip{x) , (13) 

where e'^^^ (2)('/') space-like 4- vectors on the wave surface ip{x) which are independent on the group variable 

a; the symbols A, Ba and ipa are some constants in the Minkowski space-time. They are determined via the initial 
condition of the studied problem. It is obvious that the function (p{x) can be taken so that the field Af,{x) will be 
localized in the confined region of space. 

III. FERMIONS IN THE EXTERNAL YM FIELD 

To obtain the fermion field ^(a;) we go from Eq.(4) to the so-called quadric Dirac equation which has the following 
form: 

-m'+ {rAlTaY + 2tg (t^A^TJ (7^^^) + tg (Yk.) (^A'^.Ta) } ^x) = 0; 

f (da ~ iq ■ Af,(x)Ta) + m ] 
^{x) = I ^ " ' ^' I '^{x) (14) 

First, to derive the solution of the last equation we simplify the third term in the left-hand side of Eq.(14). 
Let the initial conditions be so that the phases in Eq.(13) satisfy the equations: 

d^-" cos((p, - (^b) = 0. (15) 
Then, using Eqs.(6)-(8) and relations for the 7-matrixes[2^ H^] we obtain after direct calculations: 
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a^-^ = i (7^"^- - 7-7^) (16) 

We note, that the term containing (fabcT'^a^'^AI^A^) disappears since the vectors e^^-j ^2); which are in the expansion 
(13), are independent on the group variable a. The same takes place with respect to all terms containing cr^'' in 
Eq.(15). 

Let us find the solution of Eq.(14) in the following form: 



<^>{x) = ^,^^{x,p) = e-'P- -F^^^i^). (17) 

where F^^aif) is some multicomponent function which is the generalized Dirac spinor. It depends on both the spin 
variable a and the variable a which specifies the state of a fermion in the space of the fundamental representation of 
the SU{N) group, thereat a = 1 ^ A''; p"^ = is some 4-vector. 

We substitute ^a,a{x,p) given by Eq.(17) into the formula (14). Using relations for the 7- matrixes psl [23 |. 
independence of e^-^-j ^^2) on the group variable a in the local frame (see Eq.(13)) as well as Eq.(lO), (11), (16), we 
obtain: 



[p' - - ^'^^2jV^^^' - 2.9 (T,AX) - *5 irk.) [^TaAl) | F,,^{ip) + t (pk) f.,„(^) = 0; 

kf, = df,ip{x); (pk) = p^fc^; (18) 

where the dot over F^_a{ip) means derivative with respect to the variable ip. 

In obtaining the last equation we neglect {d^k^^l as compared with \{pk)\ (see Eq.(18)). This means that the wave 
length Xym of the YM field is unchangeable on the scale which is of the order of the de Broglie wave length of a 
fermion Xp: 



dX 



YM 



A 



YM 



Xf 



KM) I ^ 



dX 



YM 



dx 



« ^ « 1 



(19) 



The last inequality corresponds to the so called quasi-classical approximation and means that d^^k^^ = 0. The condition 
^^k^^ — can be treated as the scale invariance of the wave surface of the YM field. Thereat, the form of the wave 
surface is determined by the harmonic functions satisfying D'Lambert equation, d^d'^''if{x) — . 

In the quasi-classical approximation governed by the inequality (19) we neglect the unexplicit dependence of k^^ on 
the eikonal (p(x) and assume that only varies along the wave surface. Then, the solution of Eq.(18) can be written 
in the form: 



F^A^) = exp ( -^g' ^^^^^^ j exp <j -^g 

p2 = m2; d^,k^ ^d^d''y^ix)^0; 

where u^ip) and Va are spinors which are elements in the space of the corresponding representations. 

The second exponent in Eq.(20) is the operator acting on the spinors Ua{p) and Va- The transforming the exponent 
in the same way as it has been done in Ref.flEl (see also Appendix A), we obtain: 
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exp I ~ig Ta- 



(pk) 



} = 



g 1 I . tan 9 ^ g 9 — tan 9 



2{pk) 



tan( 



T 



(pk) 2N 



9 (pk) 



n J V (Ay) J dcp' {Ayn 



I — f *' 



(21) 



We substitute the exponent given by the last formula into Eq.(20). After that, using Eqs.(17) and (14) we obtain 



(N'^-1)A'^ 

*<T,a(a;,j3) = <^a,a{X'P) = COS 6* • exp ( -ig'^ 2N{pk) ^ ~ 



1 - igTa 



tan 6* 
9{pk) 



>(7'^fc.)(7^Ag) 
2(pk) 



tan 9 g 1 I . tan 9 g 9 — tan ( 

~9~ "^■p) 27V \^~9~^'{pk) W 



■n J V (Ay) J dcp' {Atpn 
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{pk) V 2N 



— / d^' {a;{^') p'^) / V (^::(^") P,) ; (5. r) = (a. a'^)v5(x) = o. 



(22) 



A. General Solution of Dirac Equation in External Field 

In order to derive the general solution of the Dirac equation we need to specify the physical sense of the spinors 
Ua{p),Va and the 4-vector p which are in Eq.(22). 

First, we require that the wave function (22) coincides with the solution of the free Dirac equation at 95 = 0. This 
means that the spinors Ua{p) satisfy the relations: 



a<"'kf,A^{<p = Q)u^{p)=Q\ u^{p)ux{p') = ±2m 5„x 5pp>\ p'^ = m^, (23) 

where u^{p) are the bispinors of the free Dirac field. Thereat, the first relation in Eq.(23) fixes the fundamental set 
of solutions of Eq.(14) which is determined by the parameter p. The plus and minus signs in Eq.(23) correspond to 
the Dirac scalar production of the spinors Ua-{p) and U(^{—p), respectively. 

Let us clarify the physical sense of the 4-vector p^ arising in Eq.(17). To do it we consider the projections of the 

momentum operator, p^ = —id'^, along the fc^-dircction and on the plane which is perpendicular to the vector k^ in 
the local frame govern by Eq.(13). For definiteness sake, we take fc'' = (1, 0, 0, 1); e^^^ = (0, 1, 0, 0); 61*2) = (0, 0, 1, 0). 

Then, the function (22) is the eigenfunction of the operators p^, p^and p° — p^ which eigenvalues are p^ and p"^ and 
p^ — p^, respectively. Thereat, the operators p^, p^and p° — p^ commutate with the Hamiltonian H: 



n = 7°7 gT^A,) + 70m + gA°; A^ = (A^; A), (24) 

This means that the the combinations = p^; P"^ = p^; = p° — of the components of the vector pi^ = {p°;p) 
introduced by Eq.(17) are the quantum numbers of the solutions of the Dirac equation (14). 
As for the spinor Va, we determine it by the relations: 
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vi vp = Sa^p; Tr{Ta) = 0; Tr{Ta n) = ^5ab (25) 



Then, the function (22) can be normahzed by the (5-function as follows: 



d^xnjx,p')'f„^^{x,p) = {27Tf6^p-p'). (26) 

Direct calculations show that ^cr,a{x,p) and $_(j^q,(x, —p) are orthogonal. In this way, it is obvious, that ^a-,aix,p) 
is the so-called positively frequency function whereas ^^cr,a{x, —p) is negatively frequency oneflGllTTj. This fact allows 
us to construct the general solution of the Dirac equation which describes the states both particles and anti-particles. 
Since the quantum number = p^; = p^; P^ = p'^ ~ p^ are linearly related to the vector p we can derive the 
general solution of the Dirac equation by combining the functions ^a-.aix,p) and $_(j,a(2;, —p)- As a result, due to 
the completeness condition (26) the general solution of Eq.(14) is 

*(a;) =^ / /Tri-L ,A aaAp)'faAx,p)+blM^-aAx,-p)} 
<T,a vV (2vr)3 1^ J 

^(^) = E / /T^-^rn ^S i^lnc^iP) ^a,aix,P) + ^Ap) *-^,a(x, -p)] , (27) 
a,a -I VV (2vr)3 1^ J 

where the symbols a\, ^{p)\h\ ^^{p) and acr^a{p)] ba.a{p) are the operators of creation and cancellation of a fermion 
(a<T,a(p);at q(p)) and anti-fermion (ba^a{p);blAP)) : respectively [H [111 . Thereat, a^AP) ^tAP>'' b^AP) ^^^d 
b\ a{p) satisfy the standard commutative relations for the fermion operators. 

IV. SOLUTION OF YM EQUATION IN THE EIKONAL APPROXIMATION 

Let us fix the state of the fermion vacuum so that the bilinear combinations of the operators of creation and 
cancellation of fermions are diagonal. Substituting ^{x) and '^{x) given by Eq.(22), (27) into the formula (3) we 
derive the following after direct calculations: 

<Pp J -1/ tan0 



2N0 



AN^ipkYo^ yV(Ag(^-)P.) M^(^) y , 

\ / 

^^"^^ , ALL, ^ o^ta-n^ ^ + ^ ^ ta^^^ f 1 I , AU, 1^ ^ 



{At{^)p,) |v(4(^')P^)j ^ {aiAP)<'='iP) + b<rAp)biAP)) 



(28) 



where the angle brackets mean averaging over the vacuum state of fermions. 

When the fermion system is homogeneous and isotropic the integrals containing the square bracket are equal to 
zero owing to the relativistic invariance (see Appendix B). Then, we substitute the current given by Eq.(28) and the 
field governed by Eq.(13) into Eq.(20). In taking into account of Eqs.(lO) (20) (that leads to cancellation of the 
first two terms in the left-hand side of Eq.(12)), we derive: 
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'2 fab" sin i^b - fc) = fab \ fc COS {(fib - (fir) + {cOS {ifib - fr) COS ((^^ - ^a)} -J^ \ Bs] (29) 



(P)a.,a (P) + (p)) , (30) 



where 



C = fab fc {cos ((/36 - ipr) COS (^, - ifia)} < 0. (31) 

The equations (29), (30) are closed with respect to the unknown quantities A and Ba- Having been solved they 
determine both the fermion and gauge field by means of Eqs.(13), (22), (27) so that the wave surface ip{x) is governed 
by the relations (10), (19), (20). 

Note that in the case of the N = 2 (when the SU{2) gauge symmetry occurs) the convolution (31), containing 
cosines, always is positive since the structure constants /^j'^ are the completely antisymmetrical tensor of the third 
rang e^^ due to the Jacob equality [2^: 



C = fab fc {cos {ipb - fr) COS {ifis-fa)} = ^ siu^ {fs - (fa) >0- (32) 

a,fc=0 

This means that in the framework of the developed model there is no self-consistent solution of the Dirac and 
Yang-Mills equations in the case of the SU{2) gauge symmetry. When the group dimension is more then N — 2 the 
structure constants /^^'^ can not be expressed in terms of the tensor e^^. As a result, it possible to fix the differences 
between phase in the convolution C so that C < 0. 

As for the coefficients Bg they satisfy the set of linear algebraical equations. The matrix of this set is symmetrical 
and, moreover, its diagonal elements are not all equal to zero. This means that the equation for Bg has the unambiguous 
solution. 

As a result, we have the following. The problem governed by Eqs.(l)-(5) has the unique solution when > 3. 
The solutions are determined by Eqs.(13), (27), (29), (30) and correspond to the eikonal consideration when the wave 
surface of the fields are determined by the equations: 



(dM^)) ■ (d^^ix)) - 0; (9^9'') ^(x) = (33) 

It follows from Eqs. (13), (27), (29), (30) that the Yang-Mills and Dirac equations have the self-consistent solution 
when the fermion current compensates the current of the gauge field which takes place due to self-interaction of such 
field. In other words, in the the framework of the developed model there is no the YM field without fermions. In 
terms of QCD this means that quarks and gluons can not exist separately in such approach. 

We should note here that the second relation in Eq.(33) implies that the function f{x) which is the argument in 
the expansion (13) of the field is the so called harmonic function. Owing to the initial conditions it can be always 
taken such that the field will be localized in the confined region of space. In this way, the relations (10), (33) and 
(13) directly lead to the axial gauge (11). 

V. QUANTIZING THE YANG-MILLS FIELD 

We note, that the formula (13) determining the YM field can be rewritten as follows: 

A'^{x) = —^=J-='Y] {caiqje" exp{-iipa) exp{-iqx) + c*^{q)e*'' exp{iipa) exp{iqx)} + Bak'^ 
\'V\/2uj "T^ 

q,a 

^^^;^exp(i qx~iip{x))\ 
e-^e, = e*''e: = 0; = 1; 'Z = (w;'?); 9% = 0, (34) 
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where a = 1 ^2 takes into account of two polarizations of the YM field, V is the normalizing volume. 

Let us change the coefficients Ca{q) and c*(g) by the operators of cancellation Ca-{q) and creation cj.(g) of a 
quant of the YM field so that Ca{q) and cj.(g) satisfy the Bose-Einstein commutative relation. In the quasi-classical 
approximation (19) this means that Caio) and cl{q) commutate each other: 

Caiq) ciicf) « ciiq) Ca{q), (35) 

where a is the spin variable. 

Taking into account of Eqs.(lO) and (20), we substitute Eq.(34) into the formula (12). As a result, we derive the 
relation between the occupancy number of fermions and YM quants: 



(36) 



where all notations are the same as they are in Eq.(30), (31), (34). 

Eqs.(27), (35) allow us to obtain the energy E of interacting fermions. Calculating the energy-momentum tensor 
T^^ we derive 

E ^Yl^iq) {cl{q) C,{q)] + P^°Hl^ {&lM<c.ip) + b^APj^AP)} . (37) 

q,(T cr;ct; p 

It follows from the last expression that the energy of the system of particles which consists of fermions and gauge 
quants splits on two terms. However such additivity is fictitious since the occupancy numbers of the fermion field 
depend on the value of the YM field. Formally, this manifests itself via Eq.(36) at the quasi-classical level. 

VI. DIRAC AND YANG-MILLS FIELDS AS STRONG INTERACTING MATTER 

Let us consider the solution of Eq.(29)-(30) in detail, when > 3. We assume that the phases cpa are chosen so 
that the convolution C given by Eq.(31) is negative. 

A. Fermion fields 

In principle, Eqs. (13), (27), (30) allow us to calculate both the amplitude and phases of the YM field provided 
that the occupancy number of fermions is known. However, the relation (30) is some functional equation since the 
right-hand side of it depends on the required value A via the 4-momentum = {E{p);p) of a fermion in an external 
field (see Eq.(27)) which enters into the correlators of fermi-operators. Such complicate problem is simplified and can 
be solved in a very important case when the system of fermions is some equilibrium matter whose temperature is T. 
Then, the correlators in Eq.(30) are the equilibrium occupancy numbers n{E) which are equal to [2^ : 

{aiAp)-^Ap)) = {kAp)KAp)) =^ — — ^ — — /g.^,.- ,a - <Q')^ (3^) 

1 -I- exp ( - — ^ — - 1 1 + exp I - — ^ — - 1 

where /i is the chemical potential supposed to be the same for all type of fermions; Q'^ is the mean value of a kinetic 
momentum of a fermion, Q'' — p'^ — gT"-A^, in an external field; C/^ = (1,0,0,0) is the so-called hydrodynamics 
velocity [iil. 

The functions (22) allow us to drive the mean value of the kinetic momentum of a fermion Q^^. After a direct 
calculations we get 

iQ0)2^jJ'2 + ml; ml=m' + ^'^^l~^^^' (39) 
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Since Q" can be interpretable as the mean value of the energy E oi a fermion in an external field, the last equation 
means that the interaction of a fermion with an external field leads to the re-normalization of a fermion mass, in the 
mean. 

Substituting Eqs.(27), (38) into the formula (30), we obtain: 



z 

To derive the solution of Eq.(40), first, we assume that the occupancy numbers are not too large so that the 
Boltzmann distribution is applicable to fermions. In this case the chemical potential can be expressed via the density 
of particles ng _26j . Then, we obtain from Eq. (40) at T ^ m : 

^2(^) - 1^ K^{z)- T » m, (41) 

where K^{z) are the modified Bessel functions (the McDonald functions) [27j . 
The ratio Ki{z) / K2{z) monotonically increases, so that[27|: 



Ki[z)/K2{z) = z/2; z^O 

Ki{z)/K2{z) = l- z^oo. (42) 



Then, Eq.(41) has the unambiguous solution which is 



no 



2\C\T- 



= ( "0 
- V|C|T3 



1/2 



1/3 



z < 1 



z > 1 (43) 



Since the Bolzmann approximation is correct when [nt^/T^) <^ 1, the first formula in Eq.(43) can only be used for 
calculation of the amplitude A of the YM field. Then, we derive from Eq.(40): 



A = ( ) T; - « ( ) < 1; 

(44) 



g{N^ - 1) \2\C\Ty ' T \\C\Ty 



It follows from Eq.(44) that at (noT^^) < 1 the field amplitude A is such that the effective mass of a fermion is small, 
m* T, even in the presence of the external field, i.e. the fermions in the external field remain an ultrarelativistic 
particles as before. 

On the other hand, the wave packet of an ultrarelativistic E ^ m particle does not spread out as compared with 
the case of non-relativistic particles owing to the dispersion law which is E{p) ~ p a.t E ^ m . This means that 
fermion states in the external field are single-particle ones in this case. Then, the system of fermions can be considered 
as some matter consisting of individual fermi-particles so that the interaction of them with the field results in the 
renormalization of their masses, in the mean. In this way, the individuality of a particle keeps as soon as the density 
of the matter is not too large (noT^^) < l[28l|. 

With increasing the number of fermions (or with decreasing the matter temperature) the Boltzmann approach 

becomes unsuitable. In the case of no > T^, the chemical potential /i is of the order of /i ^ n^^^ ^ T. Then, 
transforming the integral in Eq.(40) according to Ref.fl^, we obtain: 

. f 2N \^ f 2N \i 1 a ^ 



g{N^-l)J ^~V.9(^^-1)/ ^' ' T 

The last formula shows that in decreasing the temperature of matter (or in increasing its density) the mean effective 
mass of a fermion rn* in the external field is enlarged, so that the fermions become non-relativistic particles. That 
leads to delocalizing the fermion states in the space due to spreading out of a wave function. Thus, in this case the 
matter constitutes some fermi-liquid consisting of fermions with the renormalized mass. 
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B. Yang-Mills fields 

The obtained formulae (29), (30) allow us to get the tensor F^^(x) of the gauge YM field. After direct calculations 
we get from Eqs.(2), (13): 



pr-id.^)^-id,^)^. (46) 



The derived tensor F^f^ enables us to obtain the strength of both a " electric" and " magnetic" field as well as the 
energy-momentum tensor T"^^. The diagonal components of T"^^ give the energy density w and Pointing vector S: 

V"" ^ (d^if) ■ (5» = (w; S). (47) 

The last formulae show that there is no energy flux of the gauge field through any surface confining the range where 
the YM field is. This means that the role of the YM field in the considered approximation is binding fermion what 
leads to the renormalization of their masses. 

VII. DEVELOPED MODEL IN CONTEXT OF QCD 

A. Effective mass of quarks 

First, we discuss applicability of the developed model to description of the strong interacting matter generated in 
collisions of heavy ions of high energies. The quasi-classicality of the model means that the occupancy number of 
particle are large. 

In the RHIC and SPS experiments the characteristic temperature T of an equilibrium quark-gluon plasma is 
T ^ 200 ~ AOOMeV. The estimations of the initial density of energy of the plasma give that the energy density 
ui ~ 10 Gev ■ while the volume of the fireball is not less than Vq ~ 10^ F^. Then the number of particles N 
inside the fireball is of the order of 

N^^>2.5- 10\ (48) 

that is in agreement with the quasi-classical approximation. 

The gas parameter riQ^^T^^ is of the order of (n^'^T^^) ^ 1.46 3.7 at such density of the matter. On the other 
hand, the mean effective mass of a quark is of the order of 



g{N^-l)J \2\C\Ty ' T \\C\Ty 

m,^^; (j^) » 1. (49) 



V|C|T= 

It follows from the last formulae that in the intermediate range of the density of matter no ~ {gT)^ the effective 
mass is proportional to the temperature of the matter that corresponds to the result of the calculation of the thermal 
mass of a quark in the hard loop approximation [29l Isoj : 

TO, ^ g T. (50) 

B. Hadronization 



Arising the collective fermi-liquid states of fcrmions, which are governed by Eqs.(40), (45), is typical in the situation 
when there are no any channel for the particles to escape the fermion system. In the case of a quark-gluon plasma 
quarks can go out of the system due the process of the hadronization. 
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We estimate the mass of the hadron generated in the result of the hadronization. We assume that the hadronization 
is the equihbrium phase transition of the first kind. Then, the chemical potentials of the quarks /iq , gluons /ig and 
hadrons fit are: 



fig = Mff = fJ-h = (51) 



Substituting /iq = into Eq.(40) we derive: 



Ar^l.5Tc; TO*- 1.78 Tc, (52) 

where Tc is the temperature of the considered phase transition. In obtaining Eqs.(52), we set 5 = 1; |C| ~ 1. 

Let us the pious n'^ are only created in the result of the phase transition. Since 7r°'s consist of 12 quarks (including 
anti-quarks) we equalize the number of hadrons and quarks (divided by 12) at the temperature Tc- As a result, we 
obtain: 

oo oo 

"^dx f x^dx 

(53) 



2V exp (v/x2 + (M^/ T,f) - 1 I exp + (1.78)2) +l' 

where Mh is the hadron mass. In deriving the last equation we take into account that the effective mass of a quark 
is TO, — 1.78 Tc according to Eqs.(52). 

Solution of Eq.(53) gives that Mh — 0.68 Tc- If the temperature of the phase transition is — 200MeF, the 
hadron mass is Mh ~ 136Afew, that corresponds to the mass of a free pion. 

VIII. CONCLUSION 

The quasi-classical model in the gauge SU{N) field theory is considered when the YM field is assumed to be in the 
form of the eikonal wave. The self-consistent solutions of the non-homogeneous YM equation and the Dirac equation 
in the external YM field is derived. It is shown that the considered problem is solvable when the dimension of the 
gauge group > 3. Thereat, the currents generated by fermions and gauge field exactly compensate each other. 

In terms of the multi particle problem, the obtained solutions correspond to the both individual and collective states 
of fermions in matter that depends strongly on the parameters of the problem such as the density of fermions and the 
temperature of matter. As for the YM field, its amplitude appears to depend strongly on the number of fermions so 
that the field does not exist without fermions. The derived gauge field has the form of a circulatory polarized wave 
(see Eqs.(13)) which energy is concentrated in the localized region of space. Thereat, interaction of the YM field with 
fermions leads, in the mean, to the re-normalization of a fermion mass so that it enlarges with increasing the YM 
field amplitude . 

The quantum theory of the considered model is developed in the quasi-classical approximation. The energy of the 
quantized fields is obtained . It shown that the energy strongly depends on the derived relation between the occupancy 
numbers of fermions and quants of the YM field. 

The relation of the developed model to the generally accepted results in QCD is considered. In the case of the hot 
homogeneous equilibrium quark-gluon plasma the re-normalization of a fermion mass leads to arising the thermal mass 
of a quark (see Eqs.(40), (49)) which strongly depends on the matter temperature. We show that in the intermediate 
range of the density and temperature of the plasma, no — g^T^, the dependence of the quark mass on the matter 
temperature and coupling constant (see Eg. (50)) corresponds to the results of it calculations which have been made 
in the hard thermal loop approximation 1291 130| early. The hadronization as the phase transition of the first kind is 
considered. In the case of the hadronization into the lightest hadrons the calculated mass of such hadrons appears to 
be of the order of the mass of a free pion. 
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Appendix A 

We expand the second exponent in Eq.(20) in the series: 



exp < 



4 



+1 {k,^) (A^(^) 7^) dcp' (A^^iip') p,)j ^ (c^ + b) (Al) 
Taking into account of Eqs.(6)-(8), (10), (11) the even and odd terms in Eq. (A.l) can be rewritten as follows: 



^ ( jdip'{A-pf^) +iij>^k^){rA-)]\ 

exp { ~ig Ta- 



V 



{pk) 



1 / g - 2 



/ even 



2! \ (pk) V 4iV 



4! \ (pk) V 4iV 

b ( g lN^~2\ sine 



6! \ {pk) V 47V 



2 \ (pk) V 4iV 



(A2) 
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exp < 



V 



ipk) 



1 



(pk) \ 3! \ (pk) V 4^ 



{c^ + b) + 



odd 



5! Upk) y AN I 
iqTaB" (sine b 



(c^ + 2c^6)-^ 



7! I {pk) V 4:N 



(c^ + 36c4) 



{pk) \ e 2 \ {pk) V AN 



N^-2\ 6'cos6»-sin( 



^3 



(A3) 



where the parameters c and b are determined in Eq.(Al). 

Summing the exponents in Eq.(A.2), (A. 3) we directly go to Eq.(21). 



Appendix B 

Due to the relativistic invariance and Eqs.(lO), (11) the integral containing the first term in the square bracket in 
Eq.(28) is equal to zero: 



E 



d^p j gk" 
2p(o)(27r)3|(pfc) ' 



A^^ - 2 tan( 



47V e 

d^p \ 1 



^2 _ 2 , 

^'"'^ ^A^^^^'^^Ey 2p(0)(27r)3] (pfc) 
r/i(^) (AJ^(^)fc^) = 0, 



tan 



(4,a(pXa(p) + ba,a{p)bi^aiP)) = 
■ (aa,a(P)«<T,a(P) + f>a,a{p)bi,aiP)) 



(Bl) 



where fi{<f) is some scalar function. 

Because of the relativistic invariance, the fact that is a SU{N) vector, and Eqs.(lO), (11) the integral containing 
the second term in the square bracket in Eq.(28) is also equal to zero: 



E 



d^p J gk" 



2p(0)(27r)3 I {pk) 



{A^,{^)p,) j d^' {Al{^') p^) 



{pk)^e^ y 



(Oa,a(P)«a,a(p) + K,cc{jp)bl,oc{P)) = 

d^p J ppp^ p^ 9 tan^ 9 + 6 — tan ( 



'^^ V 

(4,a(PXa(p) +^<7,a(p)&i,a(p)) = 

A;X(<p) j VA^(<^') y V4(<p") (fc^fc,fc^/2(<^)+e,MfcV3(¥')+a^fc,/4(v)+a,^fc^/5(v)) =0 



(B2) 





where /i(<^) are some scalar functions. 



